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Abstract. The location of the red edge of the ZZ Ceti instability strip is defined observationally as being the lowest 
temperature for which a white dwarf with a H-rich atmosphere (DA) is known to exhibit periodic brightness 
variations. Whether this cut-off in flux variations is actually due to a cessation of pulsation or merely due to 
the attenuation of any variations by the convection zone, rendering them invisible, is not clear. The latter is 
a theoretical possibility because with decreasing effective temperature, the emergent flux variations become an 
ever smaller fraction of the amplitude of the flux variations in the interior. In contrast to the flux variations, 
the visibility of the velocity variations associated with the pulsations is not thought to be similarly affected. 
Thus, models imply that were it still pulsating, a white dwarf just below the observed red edge should show 
velocity variations. In order to test this possibility, we used time-resolved spectra of three DA white dwarfs that 
do not show photometric variability, but which have derived temperatures only slightly lower than the coolest ZZ 
Ceti variables. We find that none of our three targets show significant periodic velocity variations, and set 95% 
confidence limits on amplitudes of 3.0, 5.2, and 8.8kms _1 . Thus, for two out of our three objects, we can rule 
out velocity variations as large as 5.4 km s^ 1 observed for the strongest mode in the cool white dwarf pulsator ZZ 
Psc. In order to verify our procedures, we also examined similar data of a known ZZ Ceti, HL Tau 76. Applying 
external information from the light curve, we detect significant velocity variations for this object with amplitudes 
of up to 4kms _1 . Our results suggest that substantial numbers of pulsators having large velocity amplitudes do 
not exist below the observed photometric red edge and that the latter probably reflects a real termination of 
pulsations. 
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1. Introduction 

In spite of advances in theoretical formulations and obser- 
vational capabilities, the details of the inner workings of 
the cool H-rich pulsating white dwarfs (DAVs or ZZ Cetis) 
continue to elude us. 
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Early attempts to explain the nature of the driving 
of the pulsations of ZZ Cetis postulated that driving oc- 
curred via the k- mechanism - i.e. in a manner akin to 
that in Cepheids and <5-Scuti stars - (e.g. Dolez & Vauclair 
1981; Dziembowski & Koester 1981; Winget et al. 1982) 
thereby ignoring the effects of pulsation on the convec- 
tion zone (also known as the "frozen-in" approximation) 
even though radiative flux transport is negligible in these 
regions. Brickhill (1991) realised that the response of the 
convection zone to the perturbation is almost instanta- 
neous (~ Is) compared to the periods of the g- modes 
(hundreds of seconds). He found that as a result, the con- 
vection zone itself can drive the pulsations: part of the 
flux perturbations entering the convection zone from the 



2 



Kotak: Pulsating or not? 



largely radiative interior are absorbed by the convection 
zone and are released half a cycle later. This interaction 
between the perturbation and the response of the convec- 
tion zone drives the g-modes. 

Within the framework of this "convective-driving" 
mechanism, recently confirmed analytically by Goldreich 
& Wu (1999a), the blue edge of the instability strip i.e. 
the hottest temperature for which pulsations are excited 
and at which pulsations are expected to be observable, is 
set by the condition that ujt c w 1 for radial order, n = 
1 and spherical degree I — 1; u is the radian frequency 
and t c the thermal time constant of the convection zone. 
(Note that r c is different from the 'global' thermal time 
scale, tt in Brickhill (1991). Where necessary, we follow 
the notation of Goldreich & Wu 1999a). 

The location of the red edge is less clearly defined. 
As long as a mode is driven, its intrinsic amplitude is 
likely to remain roughly constant as it is most likely set 
by parametric resonance with stable daughter modes (Wu 
& Goldreich 2001). As the white dwarf cools, the depth of 
the convection zone increases, and damping in the shear 
layer at the base of the convection zone becomes stronger. 
At some stage, damping will exceed the driving due to the 
convection zone. This sets the physical red edge, beyond 
which pulsations are no longer excited. 

Observationally, though, the red edge may appear at 
higher temperatures. This is because, as the convection 
zone deepens, the amplitude of flux variations at the pho- 
tosphere becomes an ever smaller fraction of the amplitude 
in deeper layers; the convection zone acts as a low-pass, 
frequency-dependent filter leading to a reduction in flux 
given by (Goldreich & Wu 1999a) 

(—) = (—) 1 CD 

V f J ph \ f J h yi+lu^F 

where the subscripts "ph" and "b" refer to the photo- 
spheric flux variation and that at the base of the convec- 
tion zone respectively, for a particular eigenmode. From 
Eq. (1) one expects that the observed amplitude of a given 
mode will decrease smoothly with increasing r c , and thus 
with decreasing temperature. This might seem inconsis- 
tent with the observations, which show a rather sharply 
defined red edge. However, r c depends very steeply on tem- 
perature, and hence it may also be that the observations 
do not yet resolve the transition from easily to barely vis- 
ible (see Fig. 8 in Wu & Goldreich 1999). 

Given the uncertainties in the physical processes lead- 
ing to damping, it is not clear whether the observed red 
edge corresponds to the physical red edge, or whether it 
is just an apparent red edge where pulsations are still 
driven, but do not give rise to observable flux variations at 
the photosphere. The implication, then, is that a perfectly 
constant white dwarf might actually still be pulsating! 

Theoretical uncertainties aside, there are also uncer- 
tainties in interpreting the observations. Indeed, over the 
years, the location and extent of the observed ZZ Ceti 
instability strip have undergone several transformations 



(e.g., Greenstein 1982; Robinson et al. 1995; Giovannini et 
al. 1998). This is due in part to the difficulty in accurately 
determining the atmospheric parameters of objects hav- 
ing convectively unstable atmospheres, as some version of 
the mixing length prescription usually has to be assumed. 
Additionally, at optical wavelengths, the Balmer lines at- 
tain their maximum strengths in the middle of the in- 
stability strip. The unfortunate consequence is that vary- 
ing the atmospheric parameters gives rise to only slight 
changes in the appearance of the spectra, making it dif- 
ficult to uniquely fit observed spectra (e.g., Koester & 
Vauclair 1996). Supplementary constraints in the form of 
uv spectra, parallaxes, gravitational redshifts, have almost 
become a prerequisite. 

Most studies to date have only considered flux varia- 
tions. However, velocity variations are necessarily associ- 
ated with these flux variations and though small - of the 
order of a few kms -1 - have nevertheless been measured 
in ZZ Ceti white dwarfs: securely in ZZ Psc (Van Kerkwijk 
et al. 2000), and probably also in HS 0507+0434B (Kotak 
et al. 2002). That there are negligible vertical velocity gra- 
dients in the convection zone due to damping by turbu- 
lent viscosity is a central tenet of the convective-driving 
theory (Brickhill 1990; Goldreich & Wu 1999b). Indeed 
for shallow convection zones, this simplification has been 
shown to hold in the 2D hydrodynamical simulations of 
Gautschy et al. (1996). This means that the horizontal 
velocity is nearly independent of depth in the convection 
zone so that although photometric variations become diffi- 
cult to detect around the red edge of the instability strip as 
mentioned above, velocity variations pass virtually undi- 
minished through the convection zone. 

Observationally, the consequence would be that the 
putative, perfectly constant white dwarf might reveal its 
pulsating nature by velocity variations. If objects below 
the red edge follow the same trend as the known pulsators 
i.e. longer pulsation periods and higher amplitudes with 
decreasing effective temperature (Clemens 1993), then 
we expect the highest velocity amplitudes at the longer 

600 s) periods. For ZZ Psc (a.k.a. G 29-38), a well- 
known pulsator close to the red edge, Van Kerkwijk et 
al. (2000) measured a velocity amplitude of 5.4 kms -1 for 
the strongest mode at 614 s. For all our objects, which are 
cooler, we expect velocity amplitudes at least as large as 
those measured for ZZ Psc - lower values would constitute 
a non-detection. 

Armed with this testable theoretical expectation, we 
look for variations in the linc-of-sight velocities in objects 
that lie just below the photometric red edge of the in- 
stability strip. Our three targets are chosen - subject to 
visibility constraints during the scheduled run - from the 
list of Kepler et al. (1995) who find a handful of non- 
pulsating ZZ Cetis close to the red edge of the instability 
strip. Kepler et al. (1995) inferred relatively high masses 
for two of our three targets (see Table 1). This inference, 
though dependent on the assumed convective prescription, 
would imply that these objects have instability strips at 
higher temperatures. 
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Table 1. Observations 





WD 


T eff 


logs 


V 


Start 


End 


Exposure 


No. of 


Scatter 


Measurement 




number 


(K) 






U.T. 


U.T. 


Time (s) 


frames 


(km s _1 ) 


error (kms -1 ) 


HL Tau 76 (s) 


0416+272 


11440 


7.89 ± 0.05 


15.2 


12:58:55 


14:59:06 


20 


213 


9.0 


6.7 


G 1-7 (s) 


0033+013 


11214 


8.70 ± 0.06 


15.5 


08:37:49 


10:39:50 


74 


83 


9.8 


7.7 


G 67-23 (s) 


2246+223 


10770 


8.78 ± 0.03 


14.4 


07:33:57 


08:48:15 


74 


50 


3.5 


2.4 


G 126-18 


2136+229 


10550 


8.17 ± 0.05 


15.3 


04:46:06 


06:27:25 


30 


136 


1.4 


7.6 


G 126-18 (s) 










06:58:18 


08:58:03 


74 


80 


5.2 


4.0 



Notes. The effective temperature and logg for each of the 3 red edge objects are taken from Kepler et al. (1995) and are inferred from 
optical data only using the ML2, a=l prescription for describing the convective flux. For HL Tau 76, these quantities are taken from 
Bergeron et al. (1995) and are constrained using both uv and optical spectra for ML2/ct=0.6. The scatter due to wander is calculated at a 
representative wavelength of 4341 A i.e. at H7. The measurement error is the typical internal error on the line-of-sight velocity associated 
with the fits to the line profiles as described in Sect. 4. "(s)" refers to the Oct. 1997 service run data. The number of frames refers to the 
number of useful frames. 



In addition to the above, we search for line-of-sight ve- 
locity variations in the first ZZ Ceti to be discovered, HL 
Tau 76 (Landolt 1968), using exactly the same instrumen- 
tal setup and reduction techniques so that it serves as a 
comparison case. Treating it subsequently as a ZZ Ccti 
variable (Appendix A) allows us to indirectly constrain 
the spherical degree of the eigenmodes. 

2. Observations and Data Reduction 

Time-resolved spectra of three alleged non-pulsators, G 
1-7, G 67-23, and G 126-18 (Kepler et al. 1995) and one 
known ZZ Ceti- type white dwarf, HL Tau 76, were ac- 
quired using the Low Resolution Imaging Spectrometer 
(Okc et al. 1995) on the Keck II telescope during a ser- 
vice run on the nights of 18 th and 20 th October 1997. In 
addition to the above, we took time-resolved spectra of G 
126-18 during an observing run with the same instrument 
and set-up on 11 th December 1997. Details of the obser- 
vations together with other quantities are summarised in 
Table 1. An 8'.' 7 slit was used with a 600 line mm -1 grat- 
ing, the resulting dispersion being ~1.25Apix _1 . For the 
service run, the seeing was around l'.'O, yielding a reso- 
lution of 6 A, while the December 1997 data were taken 
in slightly worse seeing conditions (1'.'2), corresponding 
to a resolution of about 7 A. A problem with the Active 
Control System meant that 22 frames of G 67-23 taken 
during the beginning of the run had to be discarded and 
4 frames taken at the very end of the run also had to be 
disregarded as the position of the target in the slit seemed 
to have altered. The last frame of HL Tau 76 was not used 
as it was clear that light was being lost. For the observa- 
tions taken during the service run, we were unfortunately 
unable to obtain a sufficient number of halogen frames so 
we decided to forgo flatfielding the spectra altogether. 

Shifts in the positions of the Balmer lines arc in- 
troduced by instrumental flexure and differential atmo- 
spheric refraction. As in Van Kerkwijk et al. (2000), we 
corrected for the former using the position of the 1 
A 5577 A sky emission line, and for the latter by computing 
the wavelength-dependent shifts using the recipe of Stone 
(1996). As the spectra were taken through a wide slit, the 
positions of the Balmer lines will also depend on the exact 



position of the target in the slit. In principle, this position 
should be tied to the position of the guide star, but in 
practice we found that there was substantial random jit- 
ter, probably due to guiding errors and differential flexure 
between the guider and the slit. 

One way of accounting for the random movements in 
the slit is to in effect 'tag' the movements of the target to 
those of a reference star that can also be accommodated 
in the slit. For the Dec. 1997 observations, the slit was 
set at a position angle such that a G-type star was in the 
slit together with G 126-18. The wandering of G 126-18 
could then be corrected for using the positions of H/3 of 
the former as fiducial points. 

An estimate of the scatter due to random wander in 
the slit can be obtained by measuring the shifts along the 
slit of the spatial profiles and translating these into scat- 
ter in the dispersion direction. To do this, we simply fit 
Gaussian functions to the spatial profiles, and fit a low 
order polynomial to the resulting centroid positions (if 
no reference star was available). The standard deviation 
of the centroid positions from the smooth curve is then a 
proxy of the scatter due to wander in the slit. We compute 
this scatter at the representative wavelength of H7. The 
resulting estimates are shown in Tabic 1. As expected, 
the use of a reference object for G 126-18 dramatically 
lowers the scatter in the measured line-of-sight velocities 
(see Table 1). However, the value for G 67-23 shows that 
commensurate results can also be obtained without a lo- 
cal calibrator if the observing conditions and guiding are 
sufficiently stable. 



3. Lightcurves 

Lightcurves for all objects were constructed by dividing 
the line-free region of the continuum between ~ 5300 - 
5700 A by the average of that region. The Fourier trans- 
forms of the light curves were calculated up to the Nyquist 
frequency and are shown in Fig. 1. No clear periodicities 
are evident for G 67-23, G 1-7, and G 126-18, as expected 
based on the limits of 0.25%, 0.37%, and 0.46% respec- 
tively set by Kepler et al. (1995). 
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Fig. 1. Fourier transforms of the light (top panels) and velocity curves (bottom panels) for all objects. The marked peaks have 
been included in the Monte Carlo simulations in Sect. 4.1. We find no significant peaks in either the flux or the velocity FTs 
for G 1-7, G 126-18, and G 67-23. All FTs are shown up to the Nyquist Frequency except that of HL Tau 76 for clarity (see 
also caption to Fig. A.l). 



4. Line of sight velocities 

Velocity curves were constructed by measuring the 
Dopplcr shifts of three Balmer lines: H/3, H7, and H5. This 
was done by fitting a sum of a Gaussian and Lorentzian 
profile and approximating the continuum by a line. The 
centroids of the Gaussian and Lorentzian functions were 
forced to be the same. The wavelength intervals chosen for 
our fits were ^388 A, 300 A, and 164 A centred approxi- 
mately on H/3, H7, and H<5 respectively. We computed the 
velocity curve by simply averaging the velocities derived 
from each line. 

The measurement errors are of the order of the scatter 
due to random wander in all but one of the series of spectra 
(Tabic 1). The one exception is the series taken of G 126-18 
in Dec. 1997, where the scatter is much smaller because 
of the tie to a reference object. Our measurements are 
most precise for G 67-23. We include our estimates for the 
scatter due to random wander in the least squares fit of 



the velocity curve (see below); the resulting fits all have 
X r 2 cd - !• 

The Fourier Transforms of the velocity curves are 
shown in Fig. 1. None of the red edge objects show a clear 
peak well above the noise level at any frequency. 

For HL Tau 76, two of the stronger peaks have cor- 
responding peaks in the light curve. However, there are 
also two other peaks evident in Fig. 1 at 6651 /xHz and 
6928 ^Hz that are not present in the light curve at ampli- 
tudes greater than 0.9 ± 0.8% and 0.6 ± 0.8% respectively 
i.e. well below our detection limit. We will return to this 
point in Appendix A. 

Although it is clear from Fig. 1 that none of our objects 
show a large (ii 5 km s -1 ) peak at any frequency, given the 
disparities in the quality of our various data sets we must 
carry out tests in order to be able to place meaningful 
limits. We describe below two different Monte Carlo tests 
we devised to obtain quantitative estimates. 
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4.1. Significance of possible detections 

We first looked for periodicities in each of the velocity 
curves by concentrating on the few peaks of relatively high 
amplitude in the Fourier transform of the velocity curves 
and fitting these successively with sinusoids, the frequen- 
cies, amplitudes, and phases all being free parameters. We 
included a low order polynomial in the fit to remove slow 
variations for all objects except G 126-18 (Dec. 1997) for 
which the calibration with respect to the reference star 
will have already removed any slow variations. 

The first test we carry out aims to determine the like- 
lihood of the peaks seen in the Fourier Transforms (Fig. 
1) as being simply due to random noise. We shuffle the 
measured velocities randomly with respect to their mea- 
surement times after removing slow variations by fitting 
a 3 rd -order polynomial for all objects (except G 126-18, 
Dec. 1997). We then determine the amplitude of the high- 
est peak in the Fourier transform at any frequency be- 
tween 700 /iHz up to the Nyquist frequency. (This choice 
reflects an extended ZZ Ceti period range). We repeat this 
procedure 1000 times, counting the number of artificial 
data sets in which the highest peak exceeded our maxi- 
mum measured Ay. The higher this number, the higher 
the chance that even the strongest peak is not genuine. 
Table 2 shows our results. The strongest peaks in all the 
red edge objects can be attributed to noise. Only in HL 
Tau 76, the one known pulsator, is there an indication that 
the velocity signal is real. We show in Appendix A that 
external information (from the light curve) is necessary to 
confirm this. 

Table 2. Results from Monte Carlo simulations of velocity 
amplitudes 



Largest Peak in F.T. Detection 

Object v A v P(> A v ) Limit (~95%) 

QHz) (kms- 1 ) (%) (kms- 1 ) 

G 1-7 (s) 753 6.0 ± 1.5 42 SO 

G 126-18 6518 3.4 ± 1.0 67 5.2 

G 126-18 (s) 2430 4.1 ± 1.2 35 5.9 

G 67-23 (s) 3988 2.0 ± 0.7 60 3.0 

HL Tau 76 (s) 1933 4.1 ± 1.0 17 5.1 



Notes. A v is the maximum measured velocity amplitude for each 
object at the frequency (y) specified in the first column. P(> A v ) 
is the percentage probability that the largest peak (i.e. of amplitude 
A v ) is present solely by chance. The third column lists the velocity 
amplitude necessary to constitute a detection at the 95% confidence 
level — at any frequency — as estimated from our second Monte Carlo 
test (see also Fig. 2). Note that the frequency corresponding to the 
maximum velocity peak is different for the two G 126-18 data sets. 

Our second objective is to determine the minimum ve- 
locity amplitude in each data set for which we could have 
had a high (^ 95%) degree of confidence that, had such 
a signal been present, we would have detected it. To do 
this, we inject an artificial signal in our data of a given 
amplitude and random frequency and phase. The random 
frequencies are chosen from 700 ^Hz up to the Nyquist 
frequency, while the phases are chosen from to 27r. We 
repeat the procedure 1000 times and count the number of 
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Fig. 2. Probabilities that the velocity signals with the given 
amplitudes that would have been detected in the different data 
sets had they been present. The dashed line shows the 95% 
confidence level, at or above which measured velocities would 
be taken to be significant. For our best data set (G 67-23), 
we can exclude periodic signals with amplitudes greater than 
3 kms -1 with 95% confidence. 

artificial data sets in which the highest peak in the Fourier 
transform is located within a range corresponding to the 
input random frequency plus or minus the frequency res- 
olution of the data. We repeat the procedure, varying the 
input amplitude, until we recover the artificial peak as 
the maximum peak in at least 95% of the trials within the 
expected frequency range. We expect that the velocity am- 
plitude we obtain in this manner will be somewhat higher 
than the detection limit derived above. This is because for 
peaks having an amplitude exactly at the detection limit, 
there is an equal chance of a noise peak either enhancing 
or diminishing it. Thus for a peak exactly at the detec- 
tion limit, we cannot expect >50% confidence. This test 
implicitly assumes that none of the signal in our Fourier 
transform is real. For HL Tau 76, we show in the appendix 
that the higher peaks in the velocity Fourier Transform 
that are coincident with peaks in the Fourier Transform 
of the light curve are real. Hence, the limits for any signal 
not associated with flux variations will be slightly lower 
than those inferred from Fig. 2. 

The results of the above simulation are displayed in 
Fig. 2. We can rule out velocity amplitudes larger than 
5.4 kms -1 - as detected by Van Kerkwijk et al. (2000) for 
ZZ Psc - for all objects except G 1-7. Our best limit for a 
minimum significant amplitude (3kms _1 ) comes from G 
67-23. 

5. Discussion and Conclusions 

We have analysed velocity curves of three white dwarfs 
below the photometric red edge of the instability strip to 
check whether the observed red edge corresponds to the 
physical one. For ZZ Psc, a white dwarf close to the red 
edge but still within the instability strip, Van Kerkwijk 
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et al. (2000) found a maximum velocity amplitude of 
5.4 ±0.8 km s -1 . We had expected larger modulations but 
for two objects we can exclude the presence of signals as 
strong as those in ZZ Psc with high confidence; for one 
of these two, G 67-23, we can exclude the presence of any 
signal greater than 3kms _1 . 

The above indicates that pulsations actually cease be- 
low the observed red edge. How secure is this conclusion? 
We see three weaknesses. The first is that in principle, de- 
structive beating between different components in a mul- 
tiplet might conceal a real signal in our short data sets. 
This, however, would be unlikely if multiple real signals 
were present. Furthermore, two data sets for the same ob- 
ject (G 126-18) taken during different observing runs are 
unlikely to be adversely affected. Two of the modulations 
found in G 126-18 (viz. 220 s and 269 s) seemed to be 
consistent with being at the same frequency and of the 
same amplitude - within the errors - for the two differ- 
ent observing runs. We carried out further Monte Carlo 
simulations to determine the likelihood of finding rela- 
tively strong peaks at the same frequency in two separate 
data sets. On the basis of these simulations, we conclude 
that these peak coincidences are almost certainly due to 
chance: even for the best case, the peaks at 265 s and 269 s, 
there is a 7% probability of a chance coincidence. 

The second weakness is that our sensitivity turned out 
to be more marginal than we had hoped. We included one 
known ZZ Ceti, HL Tau 76, in the sample. Treating it as a 
red-edge object afforded a check on our methods. However, 
even for this bona-fidc pulsator, we found that the velocity 
curve on its own was not sufficient to demonstrate velocity 
variations. A hint of a real signal was present, but this 
could only be confirmed by imposing external information 
from the light curve, information we lacked for the red- 
edge objects. In the process, HL Tau 76 was added to the 
select group of ZZ Cetis for which velocity variations have 
been detected. 

The above might suggest that our experiment was sim- 
ply not sensitive enough. One has to keep in mind however, 
that for the different objects rather different sensitivities 
were reached. For one object, G 67-23, we can exclude 
the presence of velocity amplitudes as small as 3kms _1 
i.e. even smaller than the 4kms~ 1 peak corresponding to 
strongest mode in HL Tau 76. Furthermore, on theoreti- 
cal grounds alone, a velocity amplitude somewhat smaller 
than that of ZZ Psc is expected for HL Tau 76 given that it 
has a slightly higher temperature and its strongest mode, 
a slightly shorter period. This expectation is supported by 
the fact that the velocity amplitude of the dominant mode 
is in turn larger than that found for HS 0507+0434B, 
which is somewhat hotter still and has an even shorter 
dominant periodicity. On the same grounds, one would 
expect, as mentioned earlier, the red-edge objects to have 
larger velocity amplitudes than those observed in ZZ Psc, 



which, if present, would have been seen in two of the three 
objects. 1 

The third weakness is perhaps the most severe: G 1-7 
and G 67-23 have relatively high inferred masses, which 
implies that these objects would have ceased to be pulsa- 
tionally active at a higher effective temperature than that 
expected for a typical 0.6 Mq white dwarf, and are there- 
fore not that close anymore to the red edge. At the same 
time, G 126-18, which has a normal mass, is the coolest of 
our three observed targets. Note though that the derived 
effective temperatures and surface gravities depend on the 
treatment of convection. 

In summary, while not ironclad, our results indicate 
that pulsations have actually ceased below the observed 
red edge rather than having become photometrically un- 
dectable, and that the theoretical expectation of an ex- 
tended instability strip, beyond the observed red edge 
is flawed. In order to settle the issue observationally, 
more sensitive measurements on objects closer to the red 
edge would be worthwhile; theoretically, detailed hydro- 
dynamic modelling might go some way towards under- 
standing the interplay of the various physical processes in 
defining the red edge. 
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Appendix A: HL Tau 76 

The Fourier Transform of the light curve of HL Tau 76 
is characteristic of that of the ZZ Cetis (Fig. A.l). We 
fit the dominant peaks in the Fourier transform succes- 
sively with a function of the form Acos(2Trft' — <fi) + C 
where A is the amplitude, / the frequency, <p the phase, 
and C, a constant offset. In order to reduce the covari- 
ance between the amplitudes and phases, the time axis 
is computed relative to the the middle of the time scries 
(t 1 = t - 13 : 59 : 01.56 U.T.). We fixed the frequencies to 
those listed by Dolez & Kleinman (1996) in their prelimi- 
nary analysis of WET (Whole Earth Telescope; Nather et 
al. 1990) data as these should be more reliable than the 
frequencies inferred from our data set, which only covers a 
short time span. We include all potential real modes listed 
in Table 4 of Dolez & Kleinman (1996) with an amplitude 
> 0.9% except for modes at frequencies below ~ 1000 /iHz 
as this corresponds to the noisiest region of our Fourier 
transform. The relevant quantities are listed in Table A.l. 
We find that our measured amplitudes (Al) are on aver- 
age 1% higher than those from the WET analysis. 

Comparing the flux and velocity Fourier Transforms, 
we note first of all, that the highest peak in Fourier 
Transform of the velocity curve coincides with the dom- 
inant peak in the Fourier transform of the light curve. 
Other such coincident peaks also seem to be present, albeit 
at a lower amplitude. Since the two peaks at 6651 /iHz and 
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Fig. A.l. Fourier transform of the light curve (a) and (b) of 
the line-of-sight velocity variations shown up to 11900 /iHz for 
HL Tau 76. There are no peaks greater than 1% and 2.5 km s -1 
longward of 11900 /iHz for (a) and (b) respectively. The lower 
curve of the upper panel shows the residuals (offset by —0.8%) 
after fitting sinusoids at the frequencies indicated. 



6928 ^Hz in Fig. A. lb are not present in the light curve 
(at amplitudes greater than 0.9 ± 0.8% and 0.6 ± 0.8% 
respectively), we surmise that they are due to noise. 

Applying the external information available from the 
light curve to the velocity curve, we infer the line-of-sight 
velocities in HL Tau 76 by fixing the frequencies to those 
obtained from the light curve. The resulting values are 
listed in Table A.l. 

In order to subject these velocity amplitudes to similar 
tests as described in Sect. 4.1, we carry out another Monte 
Carlo simulation, this time taking into account the exter- 
nal information available from the light curve. We fit the 
shuffled data sets with frequencies fixed to those observed 
in the light curve. Repeating this process 1000 times, we 
find that the peak in the velocity FT at F3 - the strongest 
mode in the light curve - now has a 0.1% chance of being 
a random peak. Clearly, this improvement c.f. Table 2 is 
due to the imposition of additional information. 

Following the notation of Van Kerkwijk et al. (2000), 
we derive the relative velocity to flux amplitude ratios 
and phases, Rv and A$y. The strongest mode, F3, 
has a period of 541s that is intermediate between that 
of HS 0507+0434 (355 s, Jordan et al. 1998) and ZZ 
Psc (614s, Van Kerkwijk et al. 2000). The measured 
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Table A.l. Pulsation frequencies and other derived quantities from the light and velocity curves for HL Tau 76. 
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Fl 657 


1521 


3.3 ± 0.8 


-57 ± 14 


3.2 ± 1.0 


34 ± 18 


10 ± 4 


91 ± 23 


F2 597 


1675 


1.7 ± 0.8 


-91 ± 28 


0.3 ± 1.1 




2 ± 6 




F3 541 


1848 


5.6 ± 0.8 


42 ± 9 


4.0 ± 1.1 


84 ± 15 


6 ± 2 


41 ± 18 


F4 494 


2023 


2.6 ± 0.8 


70 ± 19 


3.7 ± 1.0 


47 ± 16 


11 ± 5 


-24 ± 25 


F5 383 


2614 


3.3 ± 0.8 


25 ± 14 


2.1 ± 1.0 


144 ± 27 


4 ± 2 


118 ± 31 


F6 781 


1280 


2.2 ± 0.8 


-81 ± 22 


3.1 ± 1.0 


48 ± 19 


18 ± 9 


130 ± 29 


Combinations: 
















Mode Period 


Frequency 


A L 




A v 




Rc 


A$ c 


F1+F3 297 


3369.0 


2.3 ± 0.8 


-20 ± 20 


1.4 ± 1.0 




6 ± 3 


-5 ± 26 


F3+F6 320 


3128.0 


2.0 ± 0.8 


-63 ± 24 


0.7 ± 1.0 




8 ± 5 


-24 ± 33 



Notes. Quantities derived from the light and velocity curves with frequencies fixed to those listed in Dolez & Kleinman (1996) whose 
nomenclature we follow. Al,v and &l,v are the amplitudes and phases obtained from fits to the light and velocity curves respectively. 
R v = A v / {2tt f A L ); A<I>v = $v - *i- Rc = A^ j A^A^) here n y = 2. A<J> C = *i ±J - ($1 ± *i) gives the relative phase of the 
combination mode with that of its constituent real modes. 
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Fig. A. 2. Results of the Monte Carlo simulation described in 
Appendix A i.e. fitting randomly shuffled velocity curves of HL 
Tau 76 with frequencies fixed to those obtained from the WET 
data. The triangles mark the velocity amplitudes thus derived. 
These amplitudes and their corresponding phases are listed in 
Table A.l. The figure shows that the signals at F3, F4, Fl, and 
F6 are unlikely to be due to chance. 



icant, that F6 has the largest Ry is consistent with it hav- 
ing a higher t value, as flux variations for higher I values 
suffer stronger cancellation while the line-of-sight veloci- 
ties do not (e.g. Van Kerkwijk et al. 2000). Comparison 
with the pulsation periods of ZZ Psc is justified; indeed, 
Dolez & Kleinman (1996) comment on the striking simi- 
larities in the modes observed for these two white dwarfs. 
A more definitive £-idcntifcation might be possible from 
analysing the wavelength dependence of pulsation ampli- 
tudes or from a much longer time series which is able to 
resolve multiplets within a period group. 



velocity amplitude at F3 also lies between that of HS 
0507+0434 (2.6±1.0kms _1 , Kotak et al. 2002) and ZZ 
Psc (5.4±0.8kms-\ Van Kerkwijk et al. 2000). The val- 
ues of Ry for all of the modes are consistent with each 
other and with having a spherical degree (£) of 1. However, 
we note with interest that the one mode in ZZ Psc iden- 
tified as having i = 2 by Clemens et al. (2000), was at a 
period of 776 s which is very close to the period of our F6. 
Even though the differences in Ry are not formally signif- 



